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ABSTRACT 

The detection of neutral deuterium in the low-metallicity damped Lyman-o' system at Zabs = 2.621 towards the quasar CTQ247 is 
reported. Using a high signal-to-noise and high spectral resolution (R = 60 000) spectrum from the Very Large Telescope Ultraviolet 
and Visual Echelle Spectrograph, we precisely measure the deuterium-to-oxygen ratio log N(D l)/N(0 1) = 0.74 ± 0.04, as well as the 
overall oxygen abundance, log Af(Ol)/W(Hl) = -5.29 ± 0.10 (or equivalently [O/H] = -1.99 ± 0.10 with respect to the solar value). 
Assuming uniform metallicity throughout the system, our measurement translates to (D/H) = (2.8^'?) x 10~ 5 . This ratio is consistent 
within errors (< 0.4 cr) with the primordial ratio, (D/H) p = (2.59 ± 0.15) x 10~ 5 , predicted by standard Big-Bang Nucleosynthesis 
using the WMAP7 value of the cosmological density of baryons (100 Sit,!' 2 = 2.249 ± 0.056). The D i absorption lines are observed to 
be broader than the O i absorption lines. From a consistent fit of the profiles we derive the turbulent broadening to be 5.2km s - ' and 
the temperature of the gas to be T = 8800 ± 1500 K, corresponding to a warm neutral medium. 
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1. Introduction 



cosmology: primordial nucleosynthesis - quasars: absoiption lines - quasars: individual: 



The primordial abundance of deuterium is fixed by Big-Bang 
Nucleosynthesis (BBN) and depends directly on the baryon- 
to-photon ratio 77. The (D/H) ratio hen ce provides a bary- 
ometer of choice (e.g. ISteigmanl I2007I) . Because of no (or 
only very marginal) subsequent production o f deuterium (e.g. 
Epst ein et al.ll976tlProdanovic & Fields 2003), and because this 
species is easily destroyed inside stars (astration), (D/H) is ex- 
pected to decrease monotonically with cosmic time. To probe 
this scenario, much effort has been devoted to the measurement 
of the deuterium abundance, in particular in low-metallicity en- 
vironments at high redshift (z > 2). This is done from observa- 
tions of D 1 an d H 1 absorption lines in QSO absorption systems 
dAdamsl [T976). Unfortunately, measurements are not yet possi- 
ble in the low density Ly-a forest because of the extremely low 
D 1 column densities expected (there is a five order of magni- 
tude difference in the abundance of the two hydrogen isotopes), 
so one has to turn to higher column density systems, namely 

Lyman-iimit systems (logiV(Hi) (cirr 2 ) > 17) and damped 2 . Observations and data reduction 



In this letter, we present a new measurement of deuterium 
abundance in the DLA at z a bs = 2.621 towards CTQ247 
(also called Q 04 05-443 or J 040718 -441014). This system was 
first reported by lLopez etldl (1200 lb from low resolution spec- 
troscopy observations, together with two other D LAs along the 
same line of sight (at z a bs = 2.551 and 2.595). IL"edoux et aT] 
(2003) subsequently confirmed the DLAs from high-resolution 
data but focused mainly on the z a b s = 2.595 system, in which 
they detected absorption lines from molecular hydrogen. The 
metal abundances in t he three DLAs were then studied by 
lLopez &Elli"sonl (12003b . We identified D 1 absorption lines in the 
z a bs = 2.621 system while screening the Very Large Telescope 
Ultraviolet and Visual Echelle Spectrograph (VLT/UVES) data 
archive. Interestingly, the system fulfils the optimal c riteria for 
the determination of (D/H) as proposed bv lYorkl (12002b : N(U 1) 
above 10 19 cm~ 2 , warm gas (6000-10000 K), wide separation of 
the D i-bearing component, and availability of 1 lines. 



Lyman-a systems (DLAs, logV(Hi) > 20.3). Although the lat- 
ter are known to be related to star-forming regions, the deuterium 
abundance is expected to rem ain close to the prim ordial value 
in the low metallicity regime dRomano et al.l 120061) . The main 
problem is then the small velocity separation between D 1 and 
Hi lines, Av ~ 80 kms" 1 , that causes Di lines to be easily lost 
within much stronger H 1 profiles. It may also be difficult to dis- 
cern betwe en D 1 lines and unrelated absorptions from the Ly-a 
forest (e.g.|Steigman 1994). This expl ains why few robus t mea- 
surements have been performed so far (Pet tini et al.l l2008). 



Based on archival ESO data Prgm. ID. 70.A-0017(A). 



CTQ 247 has been observed several times in the past decade us- 
ing VLT/UVES with different instrument setups and under vary- 
ing sky conditions. Here, we take advantage of UVES data ob- 
tained through programme ID 70.A-0017(A) (P.I. Petitjean) us- 
ing a 0.8" slit width and 2x2 CCD pixel binning. This homo- 
geneous dataset is composed of nine exposures totalling more 
than 12 hours of integration time under good seeing conditions 
(~ 0.8"). Dichroic #1 was used to observe with the blue and 
red spectroscopic arms simultaneously with central wavelengths 
adjusted to 390 nm and 570 nm respectively. Wavelength cali- 
bration frames were taken on target using a ThAr lamp imme- 
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diately after each science exposure. With this setup, a resolving 
power of R = 60 000 was reached. The spectra were reduced us- 
ing the UVES data reduction pipeline v4.9.5 based on the ESO 
Common Pipeline Library v5.3. 1 . Object and sky were optimally 
extracted using a properly over-sampled Gaussian or virtual pro- 
file depending o n the S /N raticQ. This extraction is based on a 
generalisation of Home's algorithm to get combined optimal sky 
and object flux estimates via^ 2 -minimisation. Pixels affected by 
cosmic ray hits and CCD defects were identified at the same 
time and removed via /c-sigma clipping. Individual spectra were 
shifted to the vacuum-heliocentric rest-frame before they were 
combined, by weighting the flux in each pixel by the inverse of 
its variance. The resulting spectrum has a S/N ratio per pixel of 
~ 55 at /l bs = 4000 A. Given the exquisite data quality, we were 
able to identify a slight inter-order background over-subtraction, 
at the level of 1% at /l ODS < 3800 A and increasing towards bluer 
wavelengths. We corrected this effet using the bottom of satu- 
rated Ly-ff forest lines as reference for the true zero level. 

3. Analysis 

We focus on the z aD . s = 2.621 DLA where we detected clear, non- 
saturated, and unblended Di absorption lines at z aDS = 2.62102, 
which defines the zero of our adopted velocity scale (see Fig . |TJ . 
This corresponds to the bluest component of the Oi profile, 
that spans ~ 180 kms -1 . Di lines associated with the other 
(v > kms ) components are in turn lost within the Hi pro- 
file. The total neutral hydrogen column density, logV(Hi) = 
20.45 + 0.10, was mea sured by fitting the damped Ly -o- and Ly- 
P lines dLedoux etalj|2006t lLopez & Elli"sonll2003h . Since the 
values from the two studies (also performed using different in- 
struments) agree within 0.02 dex, we did not attempt to redo this 
fit and adopt the above A^(H i)-value. 

We use standard multi-component Voigt-profile fitting to de- 
rive the column densities of the individual components, after 
normalising the spectrum locally around each line of interest. 
We simultaneously fitted five O i transitions (at A = 988.5, 988.6, 
988.7, 1039, and 1302) and the Siii/11526 transition over the 
whole profile, together with D 1/1937 and D 1/1949 in the bluest 
component. The Sin/11526 profile is located outside the Ly- 
a forest, non-affected by blends, and is characterised by non- 
saturated components and a high S/N ratio. We used this tran- 
sition as a reference to identify the individual components for 
the absorption lines located in the Ly-o- forest and to get a first 
guess of the redshifts and Doppler parameters. During the fitting 
process, the Doppler parameters of O i and Si n were assumed 
to be equal. This is the usual assumption in the field and corre- 
sponds to a Doppler parameter mostly dominated by turbulent 
motions. However, because deuterium is a very light species, the 
contribution of thermal broadening, b^ = y/2kT/m (where k is 
the Boltzmann constant, m the atomic mass and T the temper- 
ature of the gas), is not negligible. Therefore, the Doppler pa- 
rameter for atomic deuterium has been left free. The best fit is 
obtained for b(D i) = 10.1 kms 1 while b = 5.8 kms 1 forOi 
and Si n. This corresponds to a kinetic temperature of the gas 
T = 8800 + 1500 K (see Fig pi), as expected in Warm Neutral 
Medium (Peti tiean et alJl992tlWolfire et al.l!995l) . Temperature 
meas urements based on l ine widths have been difficult in DLAs 
(see ICarswell et al1l2012l) because of their components' blend- 
ing, and are thus generally limited to narrow lines which sam- 
ple cold material mu ch less representative of the typical DLAs 
(Petitiean ~et~ai1l2000h . 
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Fig. 1. Velocity plots of D i, O i, and Si n lines (black). The best- 
fit synthetic spectrum is superimposed in red. Vertical lines in- 
dicate the position of the different components. The profiles for 
D 1/1925 and D 1/1972 were computed using the derived param- 
eters but were not used in the fitting process, the former being 
contaminated by absorption from the Ly-o- forest and the lat- 
ter being saturated. We also note that the vertical lines and the 
velocity scale of the 1/1988 panel correspond to the strongest 
(0 1/1988.7) of the three transitions that produce the profile. 
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Fig. 2. Doppler parameter as a function of atomic mass. The 
solid curve and shaded area correspond to the best fit function 
b 2 = 2kT lm+b 2 mh and 1 <x range with b tm \, = 5.2±0.2 kms^'and 
T = 8800 ± 1500 K. 



See ESO UVES pipeline user manual vl9. 



The result of the Voigt-profile fitting is shown on Fig. Q~]and 
the corresponding column densities in Table Q] Formal errors 
given by fitting procedures (cr^) may underestimate the parame- 
ters uncertainties. Therefore, we also performed several fits after 
shifting the continuum by +0.5cr (where cr is the local flux error 
around each line) and take the difference in extreme values as 
the error due to uncertainties in the continuum placement, cr cont . 
Errors quoted in the table are estimated as <r 2 0t = cr 2 t + <r 2 ont . 

It is not possible to directly measure the H i column den- 
sity in each individual component but we ca n reasonably assum e 
that N(U i) scales directly with N(0 1) (see lTimmes et al.l[l997h . 
O i, H i an d D i are locked by favourable charge exchange reac- 
tions (e.g. iJenkins et al.ll2000h so we do not require any photo- 
ionisation correction. Our assumption therefore simply corre- 
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Table 1. Parameters of the modelled absorption profile 



z 


V 


h 




logiv ^cm ) 






(km s ) 


(km s~') 


Oi 


D i 


H i" 


2.621022 





5.8+0.1 (Oi)/ 10.3±0.6 (Di) 


14.22±0.02 


14.96±0.03 


19.52 


2.621403 


+31 


6.0±0.3 


14.59±0.04 


- 


19.89 


2.621492 


+39 


28.9+1.1 


14.71±0.03 


- 


20.01 


2.621589 


+47 


11.0+2.0 


14.05±0.16 




19.35 


2.621774 


+62 


5.3±0.5 


14.14±0.07 




19.44 


2.622204 


+98 


3.6±1.4 


13.29±0.09 




18.59 


2.622364 


+111 


5.3±3.4 


13.14+0.15 




18.44 


2.622876 


+ 153 


11.1+1.2 


13.62+0.11 




18.92 


2.623095 


+ 172 


6.0±2.1 


13.23±0.20 




18.53 


2.623196 


+ 180 


4.1±0.5 


13.12+0.25 




18.42 






Total 


15.15+0.02 




20.45±0.10 



Notes. (fl) Total MHl) from lLedoux et al.l f2006). distributed over the different components by scaling to A^O i). 



sponds to that of a uniform metallicity throughout the profil 
(here [O/H] = -1.99, with respect to the solar value from lLoddersI 
120031 Indeed, the war m phase of DLAs shows a high level o f 
chemical uniformity dProchaskal 120031: iRod rfguez et al. 2006), 
though we caution that this has not yet been tested at metallici- 
ties below l/50 th solar. During the fitting process, we therefore 
modelled the H 1 profiles by scaling the column densities with 
N(0 1). The Hi ^-values are highly degenerate, in particular in 
the central components, but consistent with those of 1 with ad- 
ditional thermal broadenings of ~ 10 4 K. The modelled H 1 Ly-y, 
Ly-5, Ly-e and Ly-8 profiles provide a very satisfactory fit to the 
data, further supporting our above assumption. 

4. Discussion 



Chemical evolution models (e.g. iRomano et ail 120061) predict 
very little astration of deuterium at the metallicity measured 
here, [O/H] w -2, so we can expect the deuterium abundance 
towards CTQ 247 to be very close to the primordial value. The 
D and O abundances are expected to be anti-correlated, because 
the former is destroyed in stars while the latter is produced dur- 
ing stellar nucleosynthesis. Here, the (D/O) ratio is much higher 
than what is measured lo cally, but also higher tha n several mea- 
surements at high-z (see He brard & Moosl 12003 and references 
therein). Therefore, this direct measurement alone already in- 
dicates verylittleastration. Furthermore, the lack of odd-even 
effect (lArnettll97ll) in this low-metallicity DL A indicates a low, 
quies cent star-formation history in the system (L6ne z~& Ellison! 
12001 . Indeed, we get (D/H) = 2.8™ x 10" 5 , which agrees 



remarkably well with the primordial value , (D/H) p = (2.59 + 
0.15) x 10" 5 , predicted by standard BBN dCoc et alj|2012l) using 



100Qi/i 2 = 2.249 from the study of cosmic mic rowave back- 
ground (CMB) anisotropies (Komatsu et al.ll20lH) . 

In Fig. [3] we compare our abundance measu rement to that in 
other high- redshift quasar ab sorpti on systems dBurles & Tvtlerl 
1998 a1lbT fCrighton et all l2004t [F umagall iet al.1 1201 it 
Kirkman et al] [20031: iLevshakov et al.l 120021: lO'Meara et alj 
200 U I2006T IPettini & Bowed I200U IPettini et al] 12008 



Srianand et al.l l2010h . The scatter between the different mea- 



surements is larger than the published errors. This is probably 
the consequence of underestimated uncertainties. Indeed, the ro- 
bustness of each (D/H) measurement has been long discussed in 



Strictly speaking, we do not even need the metallicity to be equal 



the lite rature. For example, the value obtained bv lCrighton et al] 
(2004) is dependent on the description of the gas kinematics and 
the er r or on N(H 1) is very likely underestimated (O'Meara et al. 
120061) . iPetfiniet al. (2008) also cautioned that the measurement 
towards Q 2206-199 7Pettini & Bowenll200ll) may suffer from a 
poor determination of N(D 1) because of the low S/N ratio and 
resolution achievecfl Overall, (D/H) ratios tend to be slightly 
higher than the standard BBN prediction, with a weighted mean 
<(D/H)> = 2.94 + °-^ x 10~ 5 . The same mean was obtained by 
Petti nieta D <H008), after removing values they considered not 
robust-enough and rescaling the errors to make them consistent 
with the observed dispersion. These authors subsequently used 
the (D/H) ratios to slightly correct downwards Q.t, from WMAP. 
Yet, the study of CMB anisotropies no w provides very robus t 
estimation of cosmological parameters (Koma tsu et al] l20llb . 
in agreement with b aryon acoustic oscillation measurements 
(lPercivaletal.ll2010h . It seems therefore more appropriate to 
seek an alternative explanation for the high (D/H) values, if 
the se were further co nfirmed. 

lOlive et al] (1201 2|) recently proposed cosmological evolution 
models that alter the nuclear processes during or right after BBN 
and which are able to solve the 7 Li problem while leading to a 
higher (D/H) p , possibly more consistent with the current data 
(horizontal dotted line on Fig. [3). It is then still possible that the 
astration of deuterium decreases the (D/H) ratio locally without 
necessarily being ac companied by an i ncreased oxygen produc- 
tion (as suggested bv l Fields et al.l2 001 and possibly observed by 
Srianand et al .11201 On) . Indeed, the astration efficiency depends 
on the relative stellar popula tions and the ratio of gas mass to 
the total galactic m ass (e.g. IVangioni-Flam & Audouzd 1988: 
IVangioni et al.ll201ll) . Local astration could therefore explain the 
non-primordial (D/H) ratios, even at low metallicity. 

Before further discussing the chemical evolution of deu- 
terium abundance, we note that metallicities are difficult to es- 
tablish precisely, in particular if elements other than 1 are used 
for low 7V(H 1) systems since ionisation corrections can be un- 
certain. For example, t he extremely low-metallicity claimed by 
Fumagalli et al.l d201 ll) might be underestimated because of the 
delicate assumption on the ionisation parameter. Some stellar 
processing of the gas is thus not completely excluded. No de- 
crease of the (D/H) ratio with metallicity is observed at high 



in all the components, but that the metallicity in the v : 
ponent equals the mean overall value. 



km s com- 



3 This system has been observed using HST-STIS while others were 
observed using echelle spectrographs on 8- 10 m class telescopes. 

4 Note however that the A^H 1) measurement is complex in this sys- 
tem owing to strong blending of the damped lines with other systems. 
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Fig. 3. High-redshift D/H values as a function of metallic- 
ity ([O/H] or [Si/H]). Our measurement towards CTQ 247 is 
shown as a red square. The red dotted lines (resp. vertical 
dashed lines) represent our constraint on (D/O) (resp. [O/H]). 
Other measurements from the literature are shown as filled cir- 
cles. Grey-filled ci rcles have cautionary remarks in the text 
(from left to right: ISrianand et al.l I20T0I ICrighton etaD l2004t 
Petti ni & Bowenll200lb . Green stars represent HD/2H2 measure- 
ments. The bl ue stripes repre sent the 1 <x and 2 <x prediction on 
(D/H) p from ICoc et al.1 d2012l) usin g standard BBN and Q,,/i 2 
from WMAP7 lkomatsu et aTTl201 ll) . The h orizontal do t ted lin e 
represent the primordial value proposed by lOlive et a"L1 (t2012l) . 
see text. 



redshift for the full range probed by direct (D/H) measure- 
ments up to one tenth solar (see Fig. [5), in agreement wit h 
standard chemical evolution models (e.g. iRomano et al. 2006). 
The high-metallicity end could however be probed at high-z 
in molecular-rich DLAs using [S,Zn/H] and the HD/2H2 ratio 
dNoterdaeme et al.ll2008l I20T0T) . As expected, the corresponding 
abundances are significantly lower than the primordial valueQ 
Although the difficulties due to line identification and small 
numbers of available transitions do not apply in the case of 
HD-based measurements, poorly known HD-che mistry and self- 
shield ing effects may alter these measurements (Balashev et al. 
120101) . Nevertheless, measurements perfo rmed at lower metal- 
licities (around -1.4, Ivanchi k et al.ll201()l and -1, Bal ashev et al.l 
120101) are very consistent with atomic-based measurements. The 
HD technique could then be a complementary way to measure 
the primordial value at [X/H] < — 1, while giving constraints on 
galactic chemical evolution at solar metallicities. 

5. Conclusion 

We report a new measurement of deuterium abundance in a 
low-metallicity DLA. The observed (D/H) value agrees re- 
markably well with the primordial value predicted by stan- 
dard Big-Bang Nucleosyhthesis using the baryon-to-photon ra- 



3 HD/2H2 values are however above what is expected in closed-box 
models and could then indicate signi ficant accretion of gas fro m the 
intergalactic medium, see discussion in Noterdaeme et al. (200J). 



tio derived independent ly from CMB anisotropics. As noticed 
by several authors (e . g.lPettini et a l. 2008; Ivanchik et al.l l2010: 
ISrianand et al.|[2010t lOlive et al .1120121) .the dispersion of other 
available measurements is larger than the published errors. This 
could be either observational or physical, where high values 
could possibly be explained by cosmological evolution models 
including post-BBN processing of the light elements, and where 
low values would be explained by local astration. Accurately 
measuring the abundance of deuterium at high redshift is clearly 
a difficult task and each measurement suffers from its own un- 
certainties. It is therefore important to obtain new high-redshift 
measurements -of both (D/H) and (D/O)- in systems with differ- 
ent chemical enrichment, for a better understanding of deuterium 
production during BBN and its astration due to stellar processing 
in the course of galactic chemical evolution. 
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